Relevance of Metric-Free Interactions in Flocking Phenomena 
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We show that the collective properties of self-propelled particles aligning with their "topological" 
(Voronoi) neighbors are qualitatively different from those of usual models where metric interaction 
ranges are used. This relevance of metric-free interactions, shown in a minimal setting, indicate that 
realistic models for the cohesive motion of cells, bird flocks, and fish schools may have to incorporate 
them, as suggested by recent observations. 
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Active matter, where energy is spent locally to pro- 
duce motion, such as in the cytoskeleton of living cells, 
displays a number of fascinating phenomena and raises 
important new theoretical questions 1]. One such gen- 
eral problem is that of the emergence of collective motion 
in large groups of self-propelled objects. Alongside the 
more traditional studies of animal behavior, the physics 
community focused, not surprisingly, on minimal models 
and their universal properties, following the seminal work 
of Vicsek et al. [2j]. In many of these works, models con- 
sist of interacting particles and rely on metric thresholds 
and cut-offs to delimit various interaction zones [3j. 

Recent facts, however, point towards a metric-free, 
topological determination of the "neighbors" with whom 
a given individual interacts. The StarFlag project pro- 
vided snapshots of reconstructed three-dimensional flocks 
of starlings of up to a few thousand birds. Statistical 
analysis revealed that a typical starling significantly in- 
teracts with its 7 or 8 closest neighbors, located at rather 
well-defined angular positions, but irrespective of their 
actual distance [4[. The authors then concluded to a 
"topological" determination of neighbors, although their 
results, obtained on flocks with a rather homogeneous 
density, can also be explained by a mere global rescaling 
of distances by each flock's density. Another hint at the 
relevance of metric-free interactions can be found in the 
recent work on crowd dynamics by Moussa'id et al. [5[ 
who argue that pedestrians make decisions based on the 
"angular landscape" formed by surrounding others, who 
may screen out otherwise close neighbors. Finally, such 
a screening effect is probably important in (moving) co- 
hesive groups of cells with direct cell-to-cell interactions. 

In this Letter, we investigate the influence of metric- 
free, "topological" determination of neighbors on collec- 
tive properties of self-propelled particles in a minimal 
setting. We introduce a "topological Vicsek model" in 
which constant-speed particles align their velocity with 
that of others forming their first Voronoi shell. We show 
that the introduction of Voronoi neighbors changes quali- 
tatively the collective properties of the original model. At 
finite density, true long-range order arises but, contrary 
to the metric model, no density segregation is present 



and the transition to collective motion is continuous with 
non mean-field critical exponents that we determine nu- 
merically in two dimensions. Large flocks evolving in an 
infinite domain move and spread diffusively, but quasi- 
long-range order nevertheless arises. 



In the original, "metric" , Vicsek model point parti- 
cles align locally, in competition with some noise. The 
interaction consists in taking the (polar, ferromagnetic) 
average of velocities over all particles within unit dis- 
tance. Two types of noise have been studied: the "an- 
gular" noise consists in adding a random angle to the 
orientation of the averaged velocity; with the "vectorial" 
noise, one adds a random vector to the average velocity 
before taking the orientation of the resulting vector to 
stream particles (for details, see, e.g. [7[). In all cases, 
particles perform uncorrelated random walks at strong 
noise, while they align perfectly in the absence of noise. 
An order-disorder transition happens in between, which 
was originally found continuous by Vicsek et al. [2J. If 
the existence, even in two dimensions, of a fluctuating 
collective motion phase with long-range order and alge- 
braic correlations was established early by Tu and Toner 
[6|, later studies [7H9| revealed that it is subdivided in 
two regions and that the transition is discontinuous due 
to the spontaneous emergence of dense, ordered, travel- 
ing solitary bands or sheets. The importance of these 
structures has been recognized implicitly by a number 
of groups [10( which tried to frustrate their emergence 
because they question the discontinuous character of the 
transition in the angular noise case. In the following, we 
restrict ourselves to the case of vectorial noise for which 
the discontinuous nature of the transition is admitted. 



The topological Vicsek model with vectorial noise is de- 
fined like its metric counterpart: N point particles move 
off-lattice at constant speed vq. In two dimensions (to 
which we restrict ourselves below), particle j is defined 
by its (complex) position r* and orientation 0* , updated 



at unit time steps according to 
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where £ is a random unit vector and 77 the noise strength. 

The key difference with the metric case is that the 
neighbors k are chosen to be the TV' particles forming 
the first shell around particle j in the Voronoi tessela- 
tion constructed from the particle positions at time t. 
With respect to its metric counterpart, this model has 
one less parameter since there is no interaction range. In 
most of the following, we used square domains of linear 
size L with periodic boundary conditions. Without loss 
of generality, we rescale lengthscales in order to fix the 
density p = N/L 2 = 1 . We are then left with two param- 
eters, the microscopic speed vq, which we keep constant 
at vq = 0.5 < 1/a/Jo = 1, and our main control parame- 
ter, the noise strength 77. 

At large 77, noise dominates the interaction, and one ob- 
serves essentially random configurations. Decreasing the 
noise strength, the order parameter ip(t) — |(expi#j)d 
increases to order one values (Fig.Q^). For a fixed noise 
strength in the ordered side, ip decreases slightly with 
L, but slower than a power- law (Fig. [Tb): our topolog- 
ical model, like the original Vicsek model, gives rise to 
true long-range polar order. Its ordered phase, though, 
is qualitatively different from that of the metric model: 
here, there is no coupling between local density and local 
order (Fig. [Tk), because neighbors in a sparse region are 
never disconnected, and thus low density does not nec- 
essarily induces disorder. The segregation mechanism 
which allows eventually for the emergence of bands in 
the metric model is absent. Indeed, scanning sizes up to 
L = 1024, no band-like structure was detected. 

Giant number fluctuations — a signature of fluctuat- 
ing ordered active phases [ll| — are nevertheless present, 
as in metric models: the variance An of the number of 
particles present in square boxes of linear size £ scales 
like n a with a ~ 1.75 > 1 where n — p£ 2 is the average 
number of particles in the box (Fig. QJi). This value of 
a is, however, significantly different from that known for 
the metric polar models, a = § [7j, [l4j. The exact ori- 
gin of this discrepancy, beyond the above remark about 
the absence of a basic segregation mechanism, remains 
unclear: a good reason would be the occurence of long- 
range connections, allowed by the topological determina- 
tion of neighbors, but this is not the case: distances to 
(Voronoi) neighbors are exponentially distributed in the 
ordered phase, which is qualitatively different from the 
Gaussian tail observed in the disordered phase, but both 
yield a small typical interaction range (Fig. QJ) . 

While they are no proof, the above findings, and in 
particular the absence of bands, hint that the transi- 
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FIG. 1: (color online) (a) time-averaged order parameter (p) 
vs 77 for different system sizes L. From top to bottom (the 
cyan arrow points from smaller to larger systes sizes): L = 32, 
64, 128, 256. (b) (ip) vs L at 77 = 0.6 (logarithmic scales). 
Inset: same data from which p^ — 0.741166 has been sub- 
tracted. The dashed line marks a power law decay as L~ 2 ' i . 
(c) Binder cumulant G vs noise amplitude 77 for different sys- 
tem sizes L (same symbols and colors as in (a)), (d) Giant 
density fluctuations in the ordered phase. An vs 77, at 77 = 0.6 
for different system sizes: L — 256 (circles), 512 (squares), 
1024 (blue). The dashed red line has a slope 1.75. (e) Scatter 
plot of local order parameter tpi versus local density pi com- 
puted in boxes of size £ = 16 in the ordered phase. Black dots: 
Vicsek metric model with interaction distance ro = 1 and pa- 
rameters vo — 0.5, 77 = 0.6, p — 2 and L = 256. Cyan dots: 
topological model with vo — 0.5, 77 = 0.6, p = 1 and L — 512. 
The lines passing through the middle of these two clouds of 
points are local averages in a moving window in pi. (f) (nor- 
malized) distribution P(r) of distances to Voronoi neighbours 
(L = 128). Black line: 77 = 0.8 (disordered); red line 77 = 0.55 
(ordered). Both cases yield typical distances of order r — 1. 
Inset: same in linear scales, showing that P(r) ~ r as r — > 0. 



tion of our topological model might not have the dis- 
continuous character observed in its metric counterpart. 
The behavior of the so-called Binder cumulant G = 1 — 
(ip 4 )t/(3(<p 2 ) 2 ) confirms this observation: when decreas- 
ing the noise strength 77, G exhibits a dip towards val- 
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FIG. 2: (color online) (a) Zoom on the crossing area 
of the G(rj) curves obtained at different L values (L — 
24, 32, 48, 64, 96, 128, 192, 256). (The variable log(| - G) was 
chosen to produce the rather good quality polynomial fits 
shown with solid lines), (b) A77 = \rj c (L) — rj§°| vs. L for 
77;?° = 0.61661. Circles: r) c {L\2) is the crossing point between 
the G{rf) curves of panel (a) for two consecutive sizes L\ and 
Li- Squares: n{L) is given by the position of the maximum of 
the susceptibility. The dashed line marks a power law decay 
as L~ 1,5 . (c) if x L 0,23 vs L for different 77 values using fit- 
ted order parameter curves (from top to bottom 77 = 0.6165, 
0.61658, 0.61661, 0.61664, 0.6171). The flat dashed black line 
is a guide for the eye. (d) Autocorrelation time vs. L at the 
critical point. The dashed line has a slope of 1.38. 



ues smaller than | (the expected value in the disordered 
phase) before reaching |, the ordered value (Fig. [TJ;). If 
this dip were to sharpen and deepen to negative values 
when increasing system sizes, it would be the signature 
of a discontinuous transition 12]. But here, on the con- 
trary, it gradually disappears as L is increased (Fig.[Tb), 
definitely ruling out this possibility: in the infinite-size 
limit, the Binder cumulant varies monotonically through 
the transition, which is thus continuous. To estimate the 
values of the associated critical exponents, we performed 
a classic finite-size-scaling study [13j , measuring the first 
moments of the order parameter as functions of 77 for 
various systems sizes. The correlation time t(t],L), esti- 
mated from the exponential decay of the autocorrelation 
function of the order parameter, was measured systemat- 
ically and we made sure that our runs were always longer 



than 10 t. To locate the critical point, we determined 
the 77 values at which the G{r\) curves obtained for dif- 
ferent system sizes L cross each other (Fig. [2^,). Without 
corrections to scaling, these curves should cross at the 
asymptotic critical point rj^° where the Binder cumulant 
takes a universal value G£° '. Here, the crossing points 
r/ c (Li,L2) estimated for two consecutive sizes L\ and L 2 



TABLE I: Critical 


exponent 


and Binder cumulant values. 
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0.23(3) 
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1.38(5) 


0.67(20) 


0.615(2) 



are not identical but quickly converge, as L\ 2 = \J L\L 2 
increases, to 77^° = 0.6166(1). If due to "ordinary" cor- 
rections to scaling, the convergence of the crossing points 
should be described by rj c (Li,L2) — rj^° + aL^ 2 u an d 
that of G c (Li,L 2 ) should be G C (L U L 2 ) = G~ + bL~" 
where w is the first correction to scaling exponent [13]. 
These scaling laws are well verified by our data using the 
above estimate of 77^° (Fig. [2}d, not shown for lj). Simi- 
larly, the location of the maximum of the susceptibility 
X = L d ((ip 2 ) t — (</?)() converges to 77^° with the expo- 
nent 1/V, but here from above whereas it is from below 
for the effective thresholds t] c (Li,L 2 ) determined using 
the crossings of the Binder cumulant curves. Using these 
two approaches to 77^° allows us to narrow our estimate to 
»7~ = 0.61661(3) and yields 1/V = 1.5(1) (Fig.Eb). Next, 
we verified that our finite-size estimates of the order pa- 
rameter at the critical point satisfy <£>(L,77^°) ex L~PI V 
(Fig. Ht), yielding j3/v = 0.23(3). Making use of the 
usual so-called hyperscaling relation 2/3 — 7 = vd where 
d is the space dimension, the above results allow us to 
complete our estimates of the static exponents yielding 
7/7/ = 1.54(6). Additional consistency checks of the 
above numbers (using various other quantities) were per- 
formed. For instance, the scaling of the maximum of the 
susceptibility Xmax(£) °c L 1 ^ is rather nicely verified, 
yielding an estimate of 7/1/ = 1.49(5) compatible with 
the above numbers (not shown). Finally, the divergence 
with L of the correlation time estimated at criticality, 
r(77^°, L) ex L z , yields an estimate of the dynamical expo- 
nent z = 1.38(5) (Fig. [2ji). A summary of our estimates 
is in Table H] They indicate that the critical properties of 
our model do not correspond to mean-field values (where, 
for instance /5 = i), nor to any known universality class. 
The transition to collective motion in the Vicsek topo- 
logical model and its ordered phase at finite density are 
thus different from that of its metric counterpart. Our 
model also possesses non-trivial properties in the zero- 
density limit of a large group of N particles evolving in an 
infinite domain. Even though no attractive interaction is 
present, subgroups can never detach to become indepen- 
dent, because of the topological determination of neigh- 
bors. This insures a weak cohesion: a finite group typi- 
cally expands diffusively at large times, but is neverthe- 
less able to order at low enough noise strengths (Fig. [3^) , 
something obviously impossible in metric models without 
attractive interactions. Remarkably, the order parameter 
remains statistically stationary, even on very long time 
scales, in spite of the constant expansion of the flock. 
The center of mass of the particles performs a persistent 
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FIG. 3: (color online) Ordered flock in an infinite domain (ini- 
tial conditions: random positions and random orientations in 
a unit circle), (a) mean distance (r) between Voronoi neigh- 
bors (top) and ip(i) (bottom) vs. time (TV = 2 14 , r\ — 0.61): 
after transient, the flock is ordered and spreads diffusively 
((r) ~ yt), indicated by the dashed line). The red data 
in the bottom panel is a moving average on a window of 
fixed length in ln(t). (b) Trajectory of the center of mass 
for t G [1, 10 7 ] plus snapshots of the flock in the asymptotic 
regime at t = 32 x IO 3 (red), 128 x 10 3 (green), 512 x IO 3 
(blue) 2048 x IO 3 (magenta) and 10 7 (orange) (same run as in 
(a)). The trajectory shows no sign of "aging" as the flock ex- 
pands diffusively, (c) (cp) vs TV in logarithmic scales (r) = 0.5, 
time averages over IO 7 timesteps after discarding an transient 
of IO 6 ). Dashed line: power law with exponent —0.074. 



relies on distances and involves the additional parameter 
n, so we focused on the model presented here. At fi- 
nite density, these "topological" flocking models exhibit 
collective properties qualitatively different from those of 
"metric" ones: the ordered moving phase shows novel 
long-range correlations (Fig. [TU), even though no direct 
long-range connections emerge (Fig. [if), and the transi- 
tion to collective motion exhibits critical properties dif- 
ferent from known universality classes. 

Beyond their importance in the general context of non- 
equilibrium order-disorder transitions, our results indi- 
cate the relevance of "topological" interactions at the 
collective level, following recent observations of their ex- 
istence locally. Since Voronoi neighbors remain at rel- 
atively well-defined distance from each other, our ap- 
proach could also provide a novel framework for modeling 
cohesive assemblies of cells, taking the Voronoi polygon 
around each particle as the effective extension of a "cell" . 

Nevertheless, realistic "topological" models of bird 
flocks, fish schools, or cell assemblies must incorporate 
cohesive interactions (e.g. pairwise attraction/repulsion) 
in order to account for the dynamics of large but finite 
groups in open domains. We are currently pursuing this 
line of research, together with theoretical investigations 
about the novel scaling laws uncovered here. 

We acknowledge the use of the CGAL package for con- 
structing Voronoi tesselations. Work funded in part by 
the French ANR projects Panurge and DyCoAct. 



random walk whose turns mostly reflect the fluctuations 
of the order parameter and are not accompanied by large- 
scale changes in the group structure or shape (Fig. [3p) . 
In other words, the orientation axes of the particles turn, 
but there is hardly any rotation of mass, just marginal 
changes in the shape and the Voronoi tesselation. 

Thus, if space is rescaled by a dynamical length re- 
flecting the expansion of the group (£ ~ \i), our sys- 
tem is mapped onto one where the microscopic velocity 
now goes to zero like \j\ft, evolving on a (quasi-) static 
Voronoi tessellation network: up to the rare rearrange- 
ments mentioned above, this is an XY model. Indeed, 
as for the two-dimensional XY model, only quasi-long- 
range order arises in the zero-density limit of our topo- 
logical model: Increasing N, the average order parameter 
decreases slightly, as in the finite-density case, but this 
time goes to zero as a power- law (Fig. [3}:). 

To summarize, we have introduced, motivated by re- 
cent observations on animal group motion, a "topologi- 
cal" Vicsek model in which self-propelled particles align 
with their Voronoi neighbors. We obtained similar re- 
sults for the model in which particles interact with their 
n (metrically) closest neighbors [15| . But this later case 
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